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We generalize the formalism for testing Lorentz invari- 
ance and the equivalence principle in the neutrino sector. 
While neutrino oscillation bounds constrain the region of large 
mixing of the the weak and gravitational eigenstates, we ob- 
tain new constraints on violations of Lorentz invariance and 
the equivalence principle from neutrinoless double beta decay. 
These bounds apply even in the case of no mixing and thus 
probe a totally unconstrained region in the parameter space. 



Special relativity and the equivalence principle can be 
considered as the most basic foundations of the theory 
of gravity. Many experiments already have tested these 
principles to a very high level of accuracy for ordinary 
matter - generally for quarks and leptons of the first gen- 
eration. These precision tests of local Lorentz invariance 
- violation of the equivalence principle should produce 
a similar effect - probe for any dependence of the 
(non-gravitational) laws of physics on a laboratory's po- 
sition, orientation or velocity relative to some preferred 
frame of reference, such as the frame in which the cosmic 
microwave background is isotropic. 

A typical feature of the violation of local Lorentz in- 
variance (VLI) is that different species of matter have a 
characteristical maximum attainable speed. This can be 
tested in various sectors of the standard model through 
vacuum Cerenkov radiation j^], photon decay Q, neu- 
trino oscillations p|,p|-pl| and physics 1^,0. In this 
article we extend these arguments to derive new con- 
straints from neutrinoless double beta decay. 

The equivalence principle implies that spacetime is de- 
scribed by unique operational geometry and hence uni- 
versality of the gravitational coupling for all species of 
matter. In the recent years there have been attempts to 
constrain a possible amount of violation of the equiva- 
lence principle (VEP) in the neutrino sector from neu- 
trino oscillation experiments However, these 
bounds don't apply when the gravitational and the weak 
eigenstates have small mixing. In this article we present 
a generalized formalism of the neutrino sector to test the 
VEP and point out that neutrinoless double beta decay 
also constrains the VEP. VEP implies different neutrino 
species to suffer from different gravitational potentials 
while propagating through the nucleus and hence the ef- 
fect of different eigenvalues doesn't cancel for the same 



effective momentum. Earlier results on neutrino oscilla- 
tions come out as special case from our present formal- 
ism. The main result is that neutrinoless double beta 
decay can constrain the amount of VEP even when the 
mixing angle is zero, i.e., when only the weak equivalence 
principle is violated, for which there does not exist any 
bound at present. 

We shall first present our formalism for VLI and then 
for VEP. For sake of clarity we formulate the problem 
for a two generation scenario involving v^, and ly^ with 
X = n,T,s. Neutrinos of different species may have dif- 
ferent maximum attainable velocities if there is viola- 
tion of local Lorentz invariance (VLI) and hence violation 
of special relativity [Q. We first assume that the weak 
eigenstates cannot be diagonalized simultaneously with 
the velocity eigenstates and the neutrinos are relativistic 
point particles. The effective Hamiltonian in the weak 
basis \vp vA is 



(1) 



In absence of VLI the neutrino mass matrix in the mass 
basis [vi V2] is given by 
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and the VLI part of the hamiltonian as 



H„ = 



vi 

V2 



(2) 



(3) 



to leading order in m'^ /p^. Here p denotes the momentum 
and fh the average mass, and for any quantity X we 
define SX = {Xi - X2), X = {Xi + X2)/2. 

In the absence of VLI, i.e., when the special threory 
of relativity is valid, Vi — 1, and simply becomes 
the momentum of the neutrinos. Here we are interested 
in a single neutrino beam (for neutrino oscillation experi- 
ments) or a single virtual neutrino propagating inside the 
nucleus with a particular momentum. For this reason we 
assume the momenta of both the neutrinos arc p. Then Vi 
corresponds to the maximum attainable speed of the cor- 
responding momentum eigenstates. Hence vi — V2 — Sv 
is a measure of VLI in the neutrino sector. As typical 
or "standard" maximum attainable speed 22^2^2. = ^ jg 
assumed. All previous bounds on this quantity Sv in the 
neutrino sector were derived from neutrino oscillation ex- 
periments and for that reason these bounds are valid only 
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for large gravitational mixing. As we shall point out, neu- 

trinoless double beta decay can constrain 6v even when 

the mixing angle vanishes. 

We shall not consider any CP violation, and hence 

Hm and H.^ are real symmetric matrices and Um and 

Uv are orthogonal matrices U^^ — . They can be 

X • 1 TT ( cos di sin Qi\ , „ 
parametrized as (7, = . „ /, , where w, rep- 

^ \- sm Qi cos di )' ^ 

resents weak mixing angle dm or velocity mixing angle 

d'u- We can now write down the weak Hamiltonian 

in the basis [fe in which the charged lepton mass 

matrix is diagonal and the charged current interaction is 

also diagonal, as 



1 /M+ M12Y 



Here / is the identity matrix and 

COSlQrn „ 



M± = m ± 
..2 



M12 = -- 



m 

sin 26*, 



-(5 m 



m 



J. / sin20„ 

QV 



Sm 

-— cos2(t^„ 
Am 



— sin 2(6*™ - 6'„) 
4m 



(4) 



given by M+ in (^. From (Q) it can easily be seen that 
in the region of maximal mixing cos20i, = 0, the double 
beta decay rate vanishes. Thus neutrinoless double beta 
decay doesn't constrain the amount of VLI for maximal 
mixing. However, when the mixing approaches zero, the 
most stringent bound from neutrinoless double beta de- 
cay is obtained. In this case Sv/2 can be understood as 
derivation from the standard maximum attainable speed 
V. As it is obvious, when there is no mixing the neu- 
trino oscillation experiments cannot give any bound on 
the amount of VLI, since in absence of mixing only VLI 
cannot allow neutrino oscillations. 

To give a bound on VLI in the small mixing region 
(including 9v = 9m = 0) we assume conservatively 
(m) ~ 0. We also assume dm < in, and thus |^ may 
be neglected. Due to the enhancement the nuclear 
matrix elements of the mass mechanism have to be re- 
placed by ^ ■ {Mp — Mqj,) with the nuclear radius R 
and the proton mass mp, which have been calculated in 
I p^ . Inserting the recent half life limit obtained from the 
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,25, 



Heidelberg-Moscow experiment J j^^g ^ i.;^-iU'"y, 
a bound on the amount of VLI as a function of the aver- 
age neutrino mass m can be given. The most reliable as- 
sumption for m is obtained from the cosmological bound 
m, < 40 eV 



tions, implying a bound of 

5v <4:X IQ-^^ for 



15 1, i.e., m < 13 eV for three genera- 



0. 



In the case of exact Lorentz invariance, we usu- 
ally write the mass matrix in the weak basis as 

mpp mp„, , mass mechanism of neutrinoless 



double beta decay, the decay rate 



[T°r]-i = ^Goi|A/i?P, 



(5) 



is proportional to the effective neutrino mass (m) = 
mee = Mj^. Here ME denotes the nuclear matrix ele- 
ment ME = Mp — Mgt, Gqi corresponds to the phase 
space factor defined in |l^ and m^ is the electron mass. 
The double beta observable can be written as 

< m > — U^iiTii = cos^ 0w + m2 sin^ 9^ 



However, combining the present experimental constraints 
from atmospheric and solar neutrino data as well as from 
tritium beta decay in a three neutrino framework and 
assuming a typical hierarchical mass pattern spectrum 
ma ^ mi^2 or ma ~ m2 ^ mi implies m <~ 0.08 
eV |l^ and improves the bound to 6v < 2 ■ 10^^* for 
0v = 9m = 0. 

In figure 1 the bound implied by double beta decay is 
presented for the entire range of sin^29y and compared 
with bounds obtained from neutrino oscillation experi- 
ments in [|o| . It should be stressed also that the GENIUS 
proposal of the Heidelberg group ]l7[ | could improve these 
bounds by about 1-2 orders of magnitude. 

For comparison, in the following the amount of VLI 
in neutrino oscillation experiments is calculated in this 
formalism. In the basis of the physical states I'a and Vb, 
the Hamiltonian becomes 



m -I- cos 29u 



(6) 



If mee — 0, the two physical eigenstates with eigenval- 
ues mi and m2 will contribute to the neutrinoless double 
beta decay by an amount U^imi and ?7e2™2, respectively, 
which cancels each other. However, if these two physi- 
cal states have different maximum attainable speed, cor- 
responding to VLI, this cancellation will not be exact 
for the same cut-off effective momentum in the neutrino 
propagator. As a result, even when m^e = 0, we can have 
neutrinoless double beta decay, which is proportional to 
the amount of VLI and the double beta observable is 
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where E ^ (P + ^) ^-nd 
-r-Ai? = ma — mi, 
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The new mixing angle Otot is a function of On 
the oscillation probability is now given by 



. 9 ^ „ .9 T^L 

= sm^ 29 tot sm^ — , 
A 



and 9.. and 



(9) 



where A = and Am^ = — m^. In the limit of 

vanishing neutrino masses, neutrino oscillations are im- 
plied only by VLI. In this case the oscillation probability 
becomes 

P{ve ~* ^x) = sin^ 2^1, SIT? pL5v ^ 

which corresponds to the expression obtained earlier . 
Here p denotes the total beam energy. From this expres- 
sion it beomes clear that in the case of no mixing, 9^ — Q, 
neutrino oscillation experiments don't constrain the size 
of VLI effects. 

In the following we present the formalism for violation 
of the equivalence principle (VEP). While in the final ex- 
pression the amount of VLI just will be replaced by VEP, 
the origin differs. In a linearized theory the gravitational 
part of the Lagrangian to first order in a weak gravita- 
tional field = 77^^ -I- hf_,^ (/i^^ = 2^dmg(l, 1, 1, 1)) 
can be written as £ = -^(1 + gi)/i^i/T^'', where T'"" is 
the stress-energy in the gravitational eigenbasis. In the 
presence of VEP the gi may differ. Assuming only viola- 
tion of the weak equivalence principle, the gravitational 
interaction is diagonal but the couplings differ. In this 
case there does not exist any bound on the amount of 
VEP. We point out that this region of the parameter re- 
gion is most restrictively bounded by neutrinoless double 
beta decay. 

The effective Hamiltonian in the weak basis again can 
be written as 



H=pl+ UrnHrnU„,' + UgHgU^', 

with H„i given in (@) and 
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(12) 



Compared to VLI, the expressions for this case of 
VEP remain unchanged, except for replacing Sv by ^SG. 

Again, g = Si±Bi. ^an be considered as the standard grav- 
itational couphng, for which the equivalence principle ap- 
plies. 

Thus the discussion of VLI can be directly translated 
to the VEP case and the bound from neutrinoless double 
beta decay for 9y — 9m = is now given by 

< 4 X lO^i^GeV (for to < 13eV) 
< 2 X lO^i^GeV (for m < O.OSeV) (13) 



to first order in ff? /p^. In formalisms where only vio- 
lation of the weak equivalence principle is assumed, one 
starts with Ug proportional to Uw, in which case there 
does not exist any bound from the neutrino oscillation 
experiments. 

The Hamiltonian in the weak basis is 



In this case, SG ~ ptpSg, where (j) is the background New- 
tonian gravitational potential on the surface of the earth. 
A natural choice for cj) would be the earth's gravitational 
potential (^ 10"^), but another well motivated choice 
could be the potential due to all forms of distant mat- 
ter. Unlike the case of VLI, the bound on the VEP will 
depend on what one chooses for the Newtonian potential 
(j). For this reason, here we only present the combined 
bound on (jj^g. 

In summary, we presented a general formalism for the 
study of both VLI and VEP in the neutrino sector. We 
pointed out that neutrinoless double beta decay can con- 
strain the amount of VLI or VEP. In particular, when 
the mixing of the gravitational eigenstates vanishes, the 
bounds from neutrinoless double beta decay become most 
stringent, while this region is not constrained by any 
other experiments. 
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FIG. 1. Double beta decay bound (solid line) on viola- 
tion of Lorentz invariance in the neutrino sector, excluding 
the region to the upper left. Shown is a double logarithmic 
plot in the Sv-s'm^ (26) parameter space. The bound becomes 
most stringent for the small mixing region, which has not 
been constrained from any other experiments. For compar- 
ison the bounds obtained from neutrino oscillation experi- 
ments (from [10]) in the — Vr (dashed lines) and in the 
Ve — i^fj (dashed- dotted lines) channel, excluding the region to 
the right, are shown. 
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